COSMIC SHEAR WITH THE CFHT 
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We present preliminary results of our cosmic shear survey currently in progress at the Canada- 
France-Hawaii Telescope (CFHT). We analysed 1.7 sq. degrees of high quality data (seeing 
below 1 arcsec), out of which we were able to measure a significant correlation of galaxy 
shape-orientation over several arcmin scale. We present measurements of the variance of the 
shear (7'^) and of the correlation functions {er(&)er(0)), {et{d)et{0)) , and we show that the 
signal is consistent with gravitational lensing by large scale structures predictions. The level 
of residual systematics after Point Spread Function correction is discussed and shown to be 
small compared to our measured signal. We outline several possible future evolutions of our 
work by using additional data of our lensing survey now extended to 9 sq. degrees. 



1 Introduction 

Since the pioneering works byilllillii, gravitational lensing by large scale structure has been 
recognized as a potential powerful probe of the mass distribution in the Universe at scales from 
1 arcmin to several degrees. The effect we are looking for consists of a small stretching of distant 
galaxies by the foreground matter, with an amplitude of 1-5 percents of distortion. Since the 
induced distortion is very small compared to the intrinsic ellipticities of the galaxies (~ 0.3), it 



Table 1: List of the fields. Most of the exposures were taken in the I band at CFHT. The total area is 1.7 deg^, 



and the 8 fields are uncorrelated. 



-Lcxigct 


IN cllllc 












seeing 


F14P1 


Fl 


CFH12K 




V 


5400 ser 


J-ViCXV l^UUU 


0.9" 


F14P2 


F2 


CFH12K 


764 arcmin^ 


V 


5400 sec. 


May 1999 


0.9" 


F14P3 


F3 


CFH12K 


764 arcmin^ 


V 


5400 sec. 


May 1999 


0.9" 


CFDF-03 


F4 


UH8K 


669 arcmin^ 




17000 sec. 


Dec. 1996 


0.75" 


SA57 


F5 


UH8K 


669 arcmin^ 




12000 sec. 


May 1998 


0.75" 


A1942 


F6 


UH8K 


573 arcmin^ 




10800 sec. 


May 1998 


0.75" 


F02P1 


F7 


CFH12K 


1050 arcmin^ 




9360 sec. 


Nov. 1999 


0.8" 


F02P4 


F8 


CFH12K 


1050 arcmin^ 




7200 sec. 


Nov. 1999 


0.9" 



can be observed only statistically by measuring the averaged orientation and elongation of several 
galaxies. Unfortunately telescope imaging defects (like optical distortion, tracking errors, etc..) 
induce also a significant coherent elongation of the galaxies, very similar to the gravitational 
lensing effect we wish to measure. Therefore the most challenging part of the analysis is the 
correction of the intrinsic telescope defects (which makes the Point Spread Function anisotropic) 
down to a level smaller than the expected lensi ng signal. This problem prevented the detection 
of the lensing effect in the earlier searches Recently a significant measurement of the 

galaxy alignment was found thanks to the improvement of the image quality^, but the analyzed 
fields were small, and yet it was not clear what sort of mass distribution was responsible for the 
lensing effect. 

The fundamental importaiice of a dean PSF correction pushed many groups to design specific 
methods of image analysis □'EJ'Ej'EjcII'EZI. During the preparation of our cosmic shear program, 
we tested intensively^ the Kaiser, Squires and Broadhurst (KSB) correction scheme on highly 
realistic simulations 13. We found that the KSB method can easily reach the 1% accuracy PSF 
correction even for severely deteriorated PSF's as we obtain in real observations. A 1% accuracy 
might be not sufficient to measure the gravitational lensing effect at the degree scale, but at 
the arcmin scale the gravitational lensing signal is enhanced by the non-linear evolution of the 
large-scale-structures and the variance of the shear (k^) can be as high as 3 — 5%. We are 
therefore confident that the present day technology is enough to detect the cosmic shear effect, 
even if significant improvements remain to be done to fully exploit the scientific case. 

Four independent groups almost simultaneously reported a significant detection BEISS and 
a new detection is coming with VLT datallZHl^. Here, we present our own detection c^, and 
discuss different aspects of the PSF correction and of the statistical measurement. Next we 
compare the measurements with the predictions and discuss their consistency. We show how 
our cosmic shear program will hopefully improve the situation very shortly: we outline the new 
measurements and the new tests of PSF correction we plan to do next. 

2 From the data set to the galaxy catalogues 

Table |l] show the list of the fields used in our cosmic shear detection. Despite the relatively 
large differences in the filter and the exposure time, the mean redshift depth is approximately 
1 according to the deepest spectroscopic surveys done so fari, with a dispersion probably very 
large (~ 1). It is impossible to give a more accurate determination of redshifts as we do not 
have enough colors for this. A complete scientific interpretation of our signal would require the 
missing redshift information, but in the early stages of the work we were more interested in a 
detection of the cosmic shear effect rather than its scientific exploitation. As discussed in the 
last Section, the scientific analysis requires some crucial issues to be addressed first. 
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Figure 1: The thick solid line show the measured 
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{7^), and the error bars are calculated using randomized 
ellipticity orientations. The dashed lines show the ilcr levels. From left to right: (a) the thin solid lines are 
(7^) measured in bins of A'^ galaxies, where the galaxies are chosen with respect to the star ellipticity strength ei 
and 62- The number of galaxies is then converted into a Active scale which is used to make the plot, (b) The N 
galaxies are chosen with respect to the optical distortion amplitude, (c) The A'' galaxies are chosen with respect 

to the CCD lines and columns. 



The total field covers about 6300 arcmin^ and has a number density of galaxies of about 
Ug ~ 30 gal/arcmin^. However, after a proper weighting of the galaxies their effective number 
density is about half (details are in the original papered). Our procedure of shape measurement 
usmg IMCAT0is described in details elsewherelBMil so it is unnecessary to give it here, therefore 
we assume that we have already a catalogue of PSF corrected ellipticities for the galaxies. 



3 From the galaxy catalogues to the cosmological parameters 

3.1 What is the relevant information? 

The easiest quantities to measure for the detection of the cosmic shear effect are the variance 
of the shear (7'^) and the ellipticity correlation functions. For a simplified cosmological model 
(zero cosmological constant, power law power spectrum with slope n and normalization ug, and 
a single redshift plane Zg), we can show thatu: 

(7^)V2 O.OlagJ^r^"-'' f , (1) 

V arcmin / 

where is the density of the Universe and 9 the measurement scale in arcmin. There are similar 
relations for the correlation functions, but we can define several types of correlation functions. 
Here we are interested in {et{0)et{0)), {er{6)er{0)) and {et{6)er{0)), where et and Cr are the 
tangential and the radial component of the shear respectively: 



et = -71 cos{2e gal) - 72 s'm{29gai) 

Sr = --f2 COs{2e gal) +'yiS'm{2e gal), (2) 

where ji are the Cartesian components of the shear and 9gai is the position angle of the pair 
of galaxies. The correlation functions mentioned above have a particular interest because of 
the specific signatures induced from weak lensing by large scale structures El: from the scalar 
nature of the gravity we can show that {et{0)er{0)) should vanish, that {et{6)et{0)) is positive 
and {er{6)er{0)) should become negative for a finite range of scale. 

"See Nick Kaiser's home page at |http: / /www.ifa.hawaii.edu/'-^kaiseT/ . 
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Figure 2: Measured ellipticity correlation functions (et(0)et(O)) (left) and (er(6)er(0)) (right) 



3. 2 Measurements and amplitude of the residual systematics 

We assumed that we already have a catalogue of galaxy ellipticities e corrected from the PSF 
anisotropy. We call the ellipticity dispersion. From a set of N galaxies at positions 9^ with 
ellipticities ea{Ok) we can built an estimate of the variance of the shear at the position 6i 




whose the ensemble average is {E['y'^{9i)\) = a'^ / N -\- {'j'^) . The term cr'^/N can be removed either 
by measuring it on randomized catalogues, or by choosing and unbiased estimate equal to Eq.( 
HI) without the diagonal elements. 

An estimate of the correlation function E[ei{0)ej{r)] (where i and j are either t and/or r) 
is built by summing ej(0)ej(r) for all the possible pairs of galaxies separated by a distance r 

Figure || shows the measured variance of the shear in our survey. The signal (thick line) 
exhibits the characteristic power law scale dependence as predicted from Eq.(|I]), and the am- 
plitude of the signal has the correct order of magnitude. However, as we shall see in the next 
Section, the cosmic variance and the error bars are too large to put tight constraints on the 
cosmology. 

We have shown in our original paperEl that the corrected ellipticity of the galaxies are uncor- 
related with the ellipticity of the stars, which demonstrates the low level of residual systematics 
present in the catalogues |. We have reproduced in Figure ^ (thin solid lines) the variance of the 
shear measured when the galaxies are picked up according to the local star ellipticity instead of 
taking the galaxies falling into a given smoothing window. Each estimate of the shear variance 
is done using Eq.(^, where the N galaxies have in common a similar PSF anisotropy amplitude 
on the data. It is then easy to convert to a scale 9, since the number density of galaxies 
is roughly constant. The resulting shear variance is shown as the thin solid lines on Figure 
|l] for three cases of possible source of residual systematics (see caption for details). It shows 
that any residual systematic cannot be due to star ellipticity, optical distortion and CCD frame 
alignment. 

A robust test of the gravitational lensing origin of the signal is to measure the correlation 
function as indicated above. Figure § shows the two relevant correlation functions {et{r)et{0)) 
and {er{r)er{0)). The third one, not shown here, (ej(r)er(O)) is zero as we expect from the scalar 
origin of the gravitational field. This results from the fact that systematics are almost non- 
existent in the galaxy catalogue, which is a strong supports for the cosmic origin of our signal. 

''However we found a constant bias of (ei) ~ —0.01 for all the galaxies which we corrected for. The source of 
the bias is still unclear. 
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Figure 3: The thick solid line is drawn from Figure |2| without the measurement error bars. The dashed lines show 
the prediction of (7^) from ray-tracing simulations and the cosmic variance obtained from several realizations. 
From top to bottom: {Qo,as) = (1., 1.); (0.3,0.85); (0.3,0.6). The source redshift is 1 and the power spectrum 

CDM-type. 



However the survey size is still too small to have a low noise measurement of these quantities, 
and to extract the useful cosmological information. Instead, we shall see this detection as a 
success showing the feasibility of the cosmic shear searches. 



3.3 Cosmological constraints 

Although it is hard to calculate the cosmic variance of the measured quantities at small scale, 
we can use ray tracing numerical simulations in order to coinpare our measurement with some 
realistic scenario. Such numerical simulations are available t3 and provide several realizations 
for a set of cosmological models, thus allowing a first analysis of the cosmic variance. Figure ^ is 
a comparison of our shear variance measurement (thick solid line) with three different models: 
(r^OjCTs) = (1., !•); (0.3, 0.85); (0.3, 0.6) from top to bottom (see the caption for details). The error 
bars show the cosmic variance measured out of several realizations for each model. Although we 
can only marginally reject the two extreme models, this figure shows that the data are already 
in favor of cluster normalized models, as suggested by the simplified analytical estimate Eq.(||). 



4 Future prospects 



The aim of the above Section was to present a condensed view of our recent cosmic shear 
detection with the firgt stream of data obtained at CFHT. All the relevant technical details can 
be found elsewhere here we emphasized the control of the systematics to a reasonable level 
and that the signal is consistent with the gravitational lensing predictions, although we are not 
yet in good position to put tight constraints on the cosmology. However, we now have 9 sq. 
degres in I (instead of 2), ready for scientific ajialysis, of comparable image quality of the data 
discussed here, and the numerical simulations Hj'Ej indicate that this is enough to make a very 
significant measurement of various lensing statistics. 

It is worthwhile to ask what can be done next. The next step is the confirmation of the 
cosmological nature of the signal over many different aspects. There are two areas with significant 
room for improvement: 

1-Beating the systematics: following the same strategy in our cosmic shear paper, we can 
search for systematics by binning the data in various ways. With 9 sq. degrees we can bin the 



data in a much refined way than we did before. We can study the ahgnment of the galaxies 
with the stars in bins of magnitude, size or other parameter. We can also split the survey into 
smaller parts and analyze independently the different parts. 

2-Search for specific gravitational lensing signatures: this was done with the correlation func- 
tion, but we can do much more with 9 sq. degrees. Of course the signal-to-noise of the correlation 
function will be much better, and we can study its shape as a function of size and magnitude 
of the galaxies. The variance of the shear measured with another filter than top-hat will be 
possible: for instance the Map statistic is known to have a specific angular dependence if 
induced by gravitational lensing. With a 9 sq. degrees survey we can also measure the skewness 
of the convergence which is a direct probe of the cosmological parameters n@. Such mea- 
sure would require accurate mass reconstruction in non-trivial topologies (non-straight edges 
and holes in the field), and preliminary studies show that such reconstructions give stable and 
unbiased mass maps. Peak statisticslliJ is also accessible with a good signal-to-noise with such a 
survey size. Measurement of the shear power spectrum will be done in order to check that the 
power is distributed over the different scales like what weak lensing theory predicts (and do not 
dominate at one specific scale for instance). 

It has recently been suggested that galaxies might have an intrinsic ellipticity correlationQa. 
Although this could also be tested using specific lensing statistics (yet to be developed), we can 
completely get rid of the effect by measuring the correlation between different source redshift 
planes. Since our full survey will be 16 sq. degrees in four colors, we hope to have the luxury to 
work with many different redshift planes (using photometric redshifts) and to extract without 
ambiguity the cosmic shear information and the intrinsic correlations. 
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